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Abstract 

Contracted capacity setting is a discrete and nonlinear 
optimization problem in consideration of expenses on the 
electricity from the utility and costs of the self-owned 
generating units (SOGUs) at the same time. This paper 
analyzes and compares the proposed improved Taguchi 
method and the cultured differential computation algorithm 
(CDCA) in solving this problem. The improved Taguchi 
method provides fast converging characteristics in searching 
the optimal solution through quality analysis in orthogonal 
matrices. The cultural algorithm (CA) used in the CDCA 
extracts and saves the domain knowledge or problem 
properties during the evolution process. The differential 
computation of the CDCA provides fast converging 
characteristics in searching the optimal solution through 
operations of mutation, crossover, and selection operations, 
which are efficient and different from the existing generic 
algorithm (GA). To compare the proposed methods, the 
paper employs the real data obtained from an 
optoelectronics factory in Taiwan. In comparison with the 
existing optimization methods in different optimization 
problems, the proposed approaches have superior results as 
revealed in the numerical results with respect to the 
computation time and the solution quality. 
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Introduction 

To determine the contracted capacity and the 
operation scheduling of the SOGUs is a highly 
discrete and nonlinear optimization problem. The 
discrete feature arises from the utility time-of-use 
(TOU) tariff designs for different peak, semi-peak, 
and off-peak periods of power consumption and the 
penalty for excessive load over contract, as well as the 
adjusting fees of contracted capacities. The nonlinear 
characteristics results from the operation dispatch of 



the SOGUs with nonlinear cost functions. A lot of 
approaches have been adopted in the literature to 
solve such a discrete and nonlinear optimization 
problem and some of them [l]-[9] are discussed 
below. 

In the nonlinear cost minimization problem for 
operation of a generation system, heuristic translation 
method [1] to determine the contracted capacities is 
suggested and widely used by Taiwan Power 
Company (TPC). Because of its use only in 
determination of the contracted capacity, the method 
cannot solve the nonlinear problem with the 
operations of SOGUs considered. Alternatively, to 
solve the similar optimization problems for 
operations of an energy supply system, mixed integer 
programming techniques have been used [2], [3] and 
the contracting policies for energy purchase/sale 
under liberalized energy markets [4]-[8] with 
promising results obtained in the studies. The 
framework, cost influence, optimal purchase in the 
demand bidding problem are further reviewed in [8] 
for the bidding strategy in a day-ahead market of a 
large consumer. More information can be found in [9] 
on the electricity procurement strategy in a day-ahead 
market and a subsequent adjustment market. 

The genetic algorithm (GA), the evolutionary 
programming (EP), and the Taguchi or improved 
Taguchi approaches have also been used to solve 
various optimization problems in power and energy 
systems [10]-[13], the contracted capacities setting 
problems [14], [15] for industrial customers 
with/without consideration of the SOGUs, which 
complicates the optimization problem. Required for 
the evolutionary computing for the capacity contract 
determination problem are the long searching time, 
the parents tuning effects, and the quality of the 
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solution searched, especially for the problem with 
SOGUs taken into account in this paper. The Taguchi 
or the proposed improved Taguchi approaches in [16] 
and [17], respectively, may enhance the global and 
local exploration abilities with good flexibility. 
However, the complex design for the orthogonal 
matrix included in the quality analysis may decrease 
the efficiency and the quality of outcome. 

The proposed approach of CDCA [18] may effectively 
search the best combination of contracted capacities 
and dispatch output of SOGUs. The Taguchi based 
method proposed by the authors in [17] for the 
optimal contracted capacity determination can thus 
be further ameliorated. In the proposed CDCA, the 
differential evolution (DE) [20] provides fast 
converging characteristics in searching the optimal 
solution through operations of mutation, crossover, 
and selection operations. The operations are efficient 
and different from the existing GA. The CA [18]-[21] 
in the CDCA incorporates domain knowledge or 
problem properties obtained during the evolution 
process to make the search process more efficient. 

In this paper, we have employed real data collected 
from supervisor control and data acquisition (SCADA) 
system of an optoelectronics factory with SOGUs in 
Taiwan and the practical Taiwan power system for the 
comparisons. The obtained simulation results 
demonstrate the feasibility and effectiveness of the 
proposed method in terms of convergence speed and 
solution precision compared with the heuristic 
scheme [1], GA [22], the mixed-integer linear program 
(MILP) [2], the improved Taguchi optimization [17], 
and the CDCA method [18]. 

Problem Formulation 

In this section, the current tariff structure of TPC is 
mathematically expressed with cost function of 
SOGUs (Please also refer to [17]) to clarify the discrete 
and nonlinear nature of the problem mentioned 
before. 

In the current tariff structure of TPC, the electricity 
bill of its customers is composed of the fees for 
capacity (By), energy (S/), regulated power factor (Rj), 
excessive load over the contracted capacities (Fj), and 
contracted capacity (Q). The operation costs of the 
SOGUs consist of fixed and variable costs (Q SOGUs ). 
Since the SOGUs are assumed to have been invested 
and installed, the fixed costs would not affect the 
decision of the contracted capacity and the operations 
of SOGUs, only given here are the operation costs of 



the SOGUs. On the other hand, suppose the consumer 
has the possibility of selling its excess energy (R s ). 

To avoid over-capacity penalty, auxiliary power 
sources of SOGUs are often used to supply the load 
demand exceeding the contracted capacity with extra 
operating costs. Therefore, an optimal contracted 
capacity exists with minimal total expenses of power 
consumption. The contracted capacity would 
determine the power expenses to pay the utilities and 
the related operating costs of the SOGUs. In this 
paper, by considering the expenses on the utility 
power and the operating costs of the SOGUs, the 
objective function of the optimization problem is 
formulated as in Eq.l. 

Min {J T } = J B . + S . + R t + F. + C, + C;** 71 * - R s ( 1 ) 

where Jt is the total power expenses of the consumer 
with SOGUs for the M-month planning horizon. As 
shown in Eq.l, the total power expenses Jt of the 
consumer with SOGUs include the capacity fee (By) 
and energy fee (Sy) in Eq.15 and Eq.18, respectively, 
depending on the capacities contracted with the 
utilities and the actual energies consumed by the 
consumer during different periods. Besides, as failing 
to keep the power factor above 0.8 or consuming the 
power exceeds the contracted capacity, the consumer 
has to pay for the power factor penalty (Rj) and 
excessive load penalty (Fj) thus incurred based on 
Eq.19 and Eq.24, respectively. To meet different load 
levels, adjusting the contract capacities would also 
cause the adjusting fees (Cy) in Eq.29. Plus, the 
operating costs of the SOGUs (Cy SOGUs ) in Eq.35 should 
be taken into account due to its operations to 
compensate for the power demands to avoid the 
penalty for excessive load (Fj). Without loss of 
generality, the sold excess energies generated by the 
SOGUs are considered in the last term R s of Eq.37. 

The detailed description of the items in Eq.l can be 
referred to Section 6, [17] and [1], where the discrete 
and nonlinear nature of the problem can be clearly 
observed. 

Overviews of Compared Methods 

In this section, we briefly describe referenced existing 
optimization methods, GA, PSO-integrated Taguchi 
method, and CDCA to solve the contracted capacity 
optimization problem of the power consumer with 
SOGUs. 
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The GA Method 

In the conventional GA, for genetic operations of 
selection, crossover, and mutation, usually used is the 
binary-string coding to represent a solution for 
further evolution [11]. The size of a population is, 
therefore, the number of binary strings. In the optimal 
contracted capacity problem, an integer, serving to 
represent a solution of contracted capacity, is used 
directly in a string for the operations of GA. 

1) Initialization and Selection 

For a given population size Np_ga, a set of initial Np_ga 
integer strings are created. An integer in the string is 
obtained randomly. The operator of Roulette-wheel 
selection is used. A string is selected from the current 
population to the next according to the probability 
given by: 



MC)= ,,[ M(pai,) 



(2) 



where Ugaj is the jth string of the current population, 
Nf_ga is the size of population, 
GA_g is the number of iteration, 
Fga(Ugaj) represents the fitness value of Ugaj. 

As indicated in (2), evaluated are the probabilities of 
the strings of the current population to survive in the 
next population based on their fitness. On the basis of 
the probabilities evaluated in (2), the roulette wheel 
strategy [11] is used for the random selection. 

2) Arithmetic Crossover 

In the GA, the crossover algorithm is employed as an 
arithmetic operator [11]. For each string selected 
through (2), arithmetic crossover operator is used as 
follows for each string in the population. 



TjGA_g+\ (]. v \TJ GA -g _|_ V TI GA -S 

u GA J V 1 t) u GAjselect 1 ^ / u GAjselect 2 

tjGA_s+\ j r GA_g (\-v\JJ GAj 
GAJ f u GA_seIect 1 ^ V 1 Y ) u GA_select 2 

where y is constant between [0,1]. 



(3) 



3) 



Mutation 



The strings of the new population after selection and 
crossover are evolved through the mutation operator 
with geometric mean probability [11]. It also reveals 
that the mutation probability function defined in [11] 
ensure the diversity of the searching directions as the 



population evolves to the optimal solution. The 
solution trapped in the local optimum can thus be 
avoided. 

The PSO-Integrated Taguchi Method 

As described conventional Taguchi method [16], [17], 
for each experiment, a solution of controlled variables 
(i.e. the contract capacities) is obtained for evaluation. 
Through the formulated orthogonal matrix, quality 
analysis on the particles of a population is performed 
with the costs and contribution values derived from 
the matrix for the iteration g. In the PSO-integrated 
Taguchi method [17], the particle is defined as a 
vector of the input levels of controlled variables to be 
determined as follows. 



(4) 



1 g Lfi,g > • • •> V n ,g ' • • •' FN ,g J 

where i = 1, q and n = l, ...,N, 

N is the number of controlled variables, 

q is the number of input levels as defined in the 
orthogonal matrix [16]. 

As noted, q is also the number of particles employed 
in the PSO algorithm for a population in a generation 
or iteration. Given the population of particles, the 
orthogonal matrix L U (^ N ) S can thus be formulated 

according the creation rules[16], [17]. It is noted that 
an orthogonal matrix L u (q") s corresponds to a 

generation of PSO in searching the optimal solution 
based on the contribution values with integrating the 
PSO method. 

The CDCA Method 

The cultural algorithms are evolutionary methods 
that extract information of the domain of the 
problem during the evolutionary process itself. A 
cultural algorithm contains two main parts: the 
population space, and the belief space. The 
population space consists of a set of possible 
solutions to the problem, and can be modeled 
using any population based technique, e.g. genetic 
algorithms. 

The DE [20] is a parallel direct search method with 
good convergence properties. Differential evolution is 
based on a mutation operator, which adds an amount 
obtained by the difference of two randomly chosen 
individuals of the current population. Through the 
arithmetic operators, such as mutation, crossover, and 
selection, the DE technique can find the optimal 
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solution from vectors of decision variables. 



Initialization of belief and 
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i 

Selection 



Acceptance function 
(select the accepted 
individuals to update the 
belief space) 




No 



FIGURE 1 FLOWCHART OF THE PROPOSED CDMA 

The proposed CDCA in [18] for the contracted 
capacity optimization uses the DE described above 
and integrates four control knowledge sources into 
the population space and the belief space of CA. Four 
control knowledge sources are employed in the belief 
space: situational, normative, topographical, and 
historical knowledge sources. The communication 
protocol is used to communicate between the belief 
and population spaces via acceptance and influence 
functions. The detail formulations are shown in [18]. 

The solution structure of the proposed CDCA is 
delineated in Figure 1. With the belief and population 
spaces initialized, the influence function selects the 
mutation operations based on the four knowledge 
sources Among the new population generated, the 
acceptance function then selects the individuals from 
the population space. The knowledge thus acquired is 
then used to update the belief space by determining 
the number of accepted individuals and the 
probability of the four knowledge sources to be 
applied to the next DE process. The evolutionary 
process terminates according to the given maximal 
number of generations. 

Numerical Results 

To compare and verify the effectiveness of the existing 
and proposed methods, the methods are numerically 
tested the contract capacity determination problems. 

Description of the Test Systems 

In this paper, the real data obtained from the SCADA 
system of a large optoelectronics factory with SOGUs 



in September, 2004 to March, 2005 were used. Based 
on their experiences, the load demand in this example 
can be regarded as deterministic. The case studies 
serve to investigate the effectiveness of the proposed 
method in the optimization of the contract capacities. 
Besides, since the consumer did not have the contract 
with the utilities as an IPP (independent power 
producer) to sell the energy, it was supposed there 
was no energy sold in the case studies. 

The data used include the operation costs of the 
SOGUs, the predicted power load demands, and the 
power bills from the utilities. Total seven SOGUs had 
been installed with maximum power output of 8MW 
for each. Through these methods, the optimal 
solution of the contracted capacities, the consumption 
of the power from utilities, and the total electrical 
power expenses are obtained. The algorithm was run 
in a personal computer of Pentium-4 with 3.0GHz 
CPU and 512MB RAM. 

Testing Results 

We analyze the solution quality and the efficiency, as 
well as the converging speed for performance 
comparison of the proposed and the existing methods 
described as follows. 

Method l.Use the heuristic translation method [1] to 
determine the contracted capacities and the 
operations of SOGUs (the solution) as 
adopted the factory currently, 

Method 2. Use the GA [22] approach to determine the 
solution, 

Method 3.Use the MILP method [2], 

Method 4. Use the improved Taguchi method [17], 

Method 5. Use the CDCA method as proposed in this 
paper. 

The results shown in Figure 2. are the comparison 
results of the total electrical power expenses for the 
total load demand of 276,939,539kWh during the 
seven-month planning horizon in different methods. 
The solution obtained from the proposed CDCA 
approach for the case of without SOGUs (Method 1) 
shows that the total power expenses are higher than 
all the solutions with SOGUs added. 

It is noted that though the fixed costs of the SOGUs 
are not shown in the paper, the costs of depreciation 
were taken into consideration in all the evaluations. 
As compared to the heuristic translation method [1] 
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adopted by the factory currently, the numerical 
results show that the GA [22], MILP [2], and the 
improved Taguchi [17] methods may achieve very 
close cost reductions of 10.27%(or US$1.5 million), 
10.20%(or US$1.49 million), and 12.18%(or), 
respectively, as given in Table 1. Through the 
proposed CDCA approach, saved is about US$2.1 



million or 14.56% of the total electrical power 
expenses for the seven months under consideration. It 
reveals that by implementing the proposed 
optimization approach in determining the contracted 
capacities and operations of the SOGUs during 
different periods, electricity expenses can be 
significantly saved. 



3,000,000 
2,500,000 
2,000,000 
1,500,000 
1,000,000 
500,000 




Cost (USD) 



I I I I 
I I I I I I 
I I I I I I 



i Heuristic translation method [1] 
I Genetic Algorithm [21] 
MILP [2] 

l Improved Taguchi method [17] 
CDCA [18] 



4 5 6 7 8 9 10 Month 

FIGURE 2 COMPARISONS OF THE CONTRACT CAPACITIES IN EACH MONTH FOR DIFFERENT APPROACHES 
TABLE 1 COMPARISONS OF THE TOTAL EXPENSES AND COST SAVINGS FOR DIFFERENT APPROACHES 



Method 


Heuristic translation 
method [ 1 ] 


Genetic Algorithm 
[22] 


MILP [2] 


Improved Taguchi 
method [17] 


CDCA [18] 


Power consumption 
from SOGUs (kWh) 


28,870,409 


23,051,099 


43,343,658 


94,275,541 


61,699,477 


Ratio (%) 


10.42 


8.32 


15.76 


34.29 


22.28% 


Total expenses (US$) 


14,658,627 


13,153,360 


13,163,367 


12,872,942 


12,524,258 


Cost Savings (%) 


0% 


10.27% 


10.20% 


12.18% 


14.56% 



Conclusion 

In this paper, the proposed optimization approaches 
have been proposed for determination of the 
contracted capacity. From the practical numerical 
results, the proposed CDCA approach, as compared 
with the heuristic method currently adopted by the 
factory, may save up to 14.56% of electrical power 
expenses in total. Also, in comparison with the 
existing optimization methods, such as the GA, MILP, 
and the improved Taguchi based optimization 
approaches, the proposed CDCA method has superior 
results as revealed in the numerical results in terms of 
the computation time needed and the quality of 
solution obtained. 

It is interesting to note from the obtained numerical 
results that the case without the SOGUs included had 
higher total power expenses as compared with the 
cases with SOGUs. The reason may arise from the 
short planning period (seven months) considered only, 
while the capital investments of the SOGUs and their 
infrastructures are depreciated into 20 years and 
levelized into each month in the current studies. 
Further, though no excess energy sold was assumed 
due to the regulation structure of the Taiwan power 



system, the proposed approach can be well applied to 
the cases with the selling of excess energy taken into 
account. 

Appendix 

As mentioned in Section 2, the electricity expenses of 
the consumer for the electricity from the utilities can 
be expressed, respectively, in the following 
subsections, which are quoted from [17]. 

(1) Capacity Fees (By): The fees for the contract capacity 
(Bj) consist of those for regular and reserve 
capacities, as given below. 

« « « (5) 

where 

fF, 1 xZ, 1 ,for(Z,,+Z,J-Z, 1 <0 

b jA = y„xx, 1+ (y. 2 ory,3)x (6) 

[ [(Z / 2 + Xj 3 ) -Xj J, otherwises 

Bj, 1 is the regular capacity tariff in j th month; Yj, 1, 
Yj, 2, and Yj, 3 are the tariffs of the peak, semi-peak, 
and off-peak regular capacities, respectively, in j th 
month; Xj, 1, Xj, 2, and X/, 3 are the corresponding 
contract capacities in jth month. The 1 st case of Eq.6 
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illustrates that if the sum of semi-peak and 
off-peak contract capacities Xy, 2 + Xy, 3 is smaller 
than the peak contract capacity, only needed to pay 
is the fees of peak contract capacity. It is also noted 
since in the current tariff structure Yj, 2 = Yj, 3, either 
Yj,2 or Yj,3 is used in the 2 nd case of Eq.6. By, 2 is the 
fees for reserve capacity, and 



B 



B .x 0.15 



(7) 



(2) Energy Fees (Sj): The energy fees (Sy) are based on 
the real power consumption in different periods. 



xZ.,+Y.,xZ. 

.(,1 5 j. 



+ Y. k xZ.,+Y. 1 xZ. i 



(8) 



where Yj, 4, Yj, 5, Yj, 6, and Yj, 7 are the tariffs of peak, 
semi-peak, off-peak, and Saturday semi-peak 
energy consumptions in j th month, which are the 
tariffs determined and announced by the utilities; 
Zj, 1, Zj, 2, Zy, 3, and Zy, 4 are the corresponding 
energy consumptions in p h month. 

(3) Regulated Power Factor Expenses (Rj): To regulate 
the power factor, the consumers whose average 
power factors are lower than 0.8 have to pay the 
penalties for low power factors. In contrast, those 
with higher power factors are rewarded. The 
formula used for the penalty and reward are stated 
as follows. 



R = 



(B j + S .) x 0.003 x (0.8 - p/,), pf j < 0.8 
(Bj + 5,) x 0.0015 x (0.8 - />/.), pf j > 0.8 



(9) 



where pfj is the average power factor of the 
consumer in j th month. 

fafor^-O.^O 



(4) Penalty for Excessive Eoad (Fj): The demands 
exceeding the contract capacities would be fined 
up to 2 to 3 times of the tariffs of the peak, the 
semi-peak and the off-peak, depending on the 
degree of load excess [1]. The loads exceeding the 
contract capacities can be expressed as follows. 



O =D -V X 



(10) 



where Oy, t and Dy, t, t = 1, 2, 3, are the loads 
exceeding contract capacities and the highest load 
demands of the peak, the semi-peak, and the 
off-peak, respectively. The corresponding penalties 
for different loading periods, Fj, t, t = 1, 2, 3, are 
described in Eqs.11-13. 



0,forO. 



<0 
forO<0 



<"lxV 



(11) 



137., x(Oj, -0.1xX yl ),forO yl >0.1xX., 



Eq.ll means that the excessive peak load would be 
fined at 2 times of tariff for the portion within 10% 
of the peak contract capacity, and at 3 times for the 
portion over 10%. In Eqs.12-13, to avoid double 
counting of penalties for excessive peak and 
semi-peak demands over contract capacities, only 
fined are the semi-peak demands over the peak 
ones, i.e. (O j2 — O jl )>0 and/or excessive off-peak 

demands over the larger value between excessive 
the peak and the excessive semi-peak demands, 
i.e. O., 3 -max(O. 1 ,O. 2 )>0 



27. 2 x (0. 2 - O, ,), for < (0, 2 - 0. ,) < O.lx X . 2 
|3K x|(0 . , O.lx.V .|.for(0 . ) M>. 1 x .V 



(12) 



p,forO j3 -max(O >l ,O j2 )<0 
p 3 =\2Y.^x[0 3 -max(0 . ,, 0. 2 )],for < 0. 3 -max(0. ,, 2 )<0.1xX. 3 

37 '; x {[0. 3 - max(0 .() )| O.lx.V .]. for O , - max(0 . .() . O.lx.V 



(13) 



= 1 



0,forX„<X._ u 

PMX^-X^lforX^^X., < max(X, r ) 
P, x [ max (X, ,) - X ,], for X , > max'(X, , ) 

i=l,...,y-l J J 1=1,..., j-1 



(15) 



0,forX., <X._ lt 

0,forX ; _ lf <X. t < max(X. ( ) 

Px[X - max (X.JLfor X > max (X. ,) 



(16) 



K 



0, for X , < X , , 

P 2 x[H/xX ]l -(X j _ ll +X J _ 2l +...+ X ] _ H J]-(K. ll +K. 2l +...+ K. H ^JorX 

P 2 x [H x max (X, , ) - (X , + X ,+...+ X ._ H , )] - (K jA , + K ,+...+ AT , ), for X , > max (V, , ) 

J i=l,—, J-1 ' '' i'=l,...,y'-l 



(17) 
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l*v - 


Ji ! 
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— 


! Month 



(b) 

FIGURE 1 SEVERAL EXAMPLES ILLUSTRATING THE ADJUSTING 
FEES OF_THE_CONTRACT RESERVE CAPACITY [17] 

Based on the descriptions above, the penalty for 
excessive loads i|s Expressed as in Eq.14. 



F.=F.,+Fi,+F., 

J .1,1 J, 2 ;,3 



(14) 



(5) Adjusting Fees of Contract Capacity (Cj): Based on 
the original contract capacities with the TPC, the 
consumers can adjust the contract capacities by 
paying adjusting fees, except free adjustment of 
reducing the contract capacity or shifting the peak 
contract capacity to the off-peak one [1]. The 
adjusting fees for the increased contract capacities 
can be divided into two parts comprising the 
equipment expansion fees Ej,i and E'jx and the 
equipment maintenance fees Kj.i for transmission 
and distribution (T&D) systems. 

Suppose Xj, t is the adjusted contract capacity from 
Xj-i, t, where t = 1, 2, and 3 refer to the periods of 
the peak, the semi-peak, and the off-peak, 
respectively. The calculations of Ej, t and E'j, i, as 
well as Kj, i are given in Eqs.15-17, where Pi and Pi 
are the unit prices of the equipment expansion and 
maintenance fees, respectively; max (Xi, t) is the 

consumer has ever had before for month z'=l, 
While Xj, t < max (Xi, ,), only would Ej.i be paid, 

otherwise both E/,i and E'j,\ should be paid for the 
portions with the adjusted contract capacity Xj, t 
smaller and larger than the previous maximal 
contract capacity . max , (Xi, t), respectively. Separate 

calculations of Ej, t and E'j, i are due to different 
fees for the portions of adjustments lower or 
higher than max (X .). 

The equipment maintenance fees Kj, t are 
calculated as shown in Eq.17 where Hj is the 
number of months for which the contract capacity 
Xj-H„ t is smaller than but next in value to the 
adjusted capacity X/,t in month;, i.e. 

Xj, t > Xj-H„t > Xj-i,t, i=l, 2, Hj. The instances in 



Figure 3(a) and (b), respectively, show the 
definitions of Kj, i in Eq.17 for cases of Xj-i, t < Xj, t 
max (X i , t ) and 



Xj, t > max (Xi,t). The shaded areas in Figure 3(a) 

and (b) indicate the ways how Kj, t, Kj-i, and Kj-2, i, 
for H j = 2 and 

t = l| are calculated. 

In_ _^<iditiori_ to _ the^ Eidjusting fees for contract 
capacities, Ej, i, E'j, i, and Kj, i respectively, the 
corresponding fees 

Ej, 2, E'j,' 2, and Kj, 2 of contract reserve capacities, 
should be paid, which are defined in Eq.18. 



E j2 =0.25x£. ] 



Ej 2 = 0.5 x E 
K . ={).\5xK 



(18) 



According to variation amount of the contract 
capacity, the adjusting fees Cj are calculated based 
on the T&D equipment expansion and 
maintenance fees of the contract capacity Cj,\ and 
the fees of the reserve contract capacity Cj,i, i.e. 



(19) 



where C/,i is composed of the equipment 
expansion and maintenance fees as described in 
Eq.20 [1]. 



c J*y,+*;,. K <0.5x/:' 
'•' [0.5 x Ej 4 + Ej , , otherwise 

Similarly, Cj,i is calculated below. 



c JKj^Ej^, 



otherwise 



(20) 



(21) 



(6) Operation Costs of the SOGUs (Q SOGUs ): The 
operation costs of the SOGUs consist of fixed and 
variable costs. Since the SOGUs are assumed to 
have been invested and installed, the fixed costs 
would not affect the decision of the contract 
capacity and the operations of SOGUs, only given 
here are the operation costs of the SOGUs. The 
operation costs of generating units comprise the 
fuel costs and the maintenance costs. As an 
example, in the case studies of this paper, the fuel 
costs are calculated by using Eqs.22-24. 



(22) 
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i 2 

a k,h ~ a X kh ' 



bx kh + c 



M u =^^,xP Gt)i x 0.9317 (Ukg)xC m 



(23) 
(24) 



where Xk, h is the hourly generating factor of the k >h 
unit at hour h; Pgk h (kW) and Pat, max (kW) are the 
corresponding hourly power output and the 
maximum power output of the unit, dkh (kg/kWh) 
is the hourly fuel consumption; a, b, and c are the 
constants associated with the characteristics of the 
generating unit. Mu($) is the fuel cost of the k th 
unit over the period considered, and Cm P is the 
fuel-oil price per liter in the market. In this paper, 
we use Cm V of 0.248 (US$/Z). The variable 
maintenance cost, Mi,k, are calculated for different 
generating units. The variable cost of the SOGUs, 
qsogus j s thus the sum of Mi, k and Mi, k. 



'=M Uk +M %k 



(25) 



The SOGUs studied in this paper are the 
heavy-oil-fired generators. The maximum output 
power of each generator is 8 MW, with zero 
minimal power output and fast ramping rate 
within 30 sec to start up and shut down. Compared 
with the generating cost, the start-up and 
shutdown costs of the SOGUs studied are 
relatively small and are neglected in this paper. 
The minimal up-time and down-time of the 
SOGUs are also negligible as compared to the load 
interval (one sample load per 15 minutes) and are 
thus not considered in the paper. The start-up and 
shutdown costs as well as the minimal up-time and 
down-time constraints of the SOGUs can be easily 
integrated in the model, if they are significant for 
the other types of SOGUs. 

Suppose the consumer has the possibility of selling 
its excess energy, the revenues obtained can be 
calculated by 



R s =E s xC s 



(26) 



where Es is the sold excess energy generated by the 
SOGUs in MWh, and C s is the unit price ($/MWh) 
for the transaction. 
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